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Abstract 

Through an extended kinetic model, we study the nonhnear generation of quasi-static magnetic 
fields by high-frequency fields in a plasma, taking into account the effects of the electron spin. It 
is found that although the largest part of the nonlinear current in a moderate density, moderate 
temperature plasma is due to the classical terms, the spin may still give a significant contribution 
to the magnetic field generation mechanism. Applications of our results are discussed. 
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I. INTRODUCTION 



The low-frequency fields nonlinearly generated by high-frequency waves in a plasma has 
been the subject of much studies during several decades il-15|. The focus has been on 
;he ponderomotive force, see e.g. jl, 2|, on heating nonlinearities due to collisional effects 
sl, or on the possibility to generate quasi-static magnetic fields j4-15|. Generation of such 
fields has proven to be rather sensitive to small deviations from the classical collisionfree 
Vlasov model. Thus Ref. jll| found that a small collision frequency (described by the 
Lorenz collision model) was sufficient to significantly alter the predictions from the Vlasov 
equation. Furthermore, Refs [isl found that weak relativistic effects (with the thermal 
energy much smaller than the electron rest mass energy), also could alter the nonlinear 
generation of quasi-static magnetic fields. 

n the present paper we will investigate to what extent quantum mechanical effects 
26 1 , and in particular electron spin effects 23-[30|, can change the nonlinear low- frequency 



23 



16 



response to high-frequency fields. A semi-classical model for spin effects generalizing the 



28 



Vlasov equation was presented decades ago 
magnetic dipole force and the magnetization current due to the spin 
quantum mechanical model was given 



Recently it was irnproved to include the 



29|, and finally a fully 



30] , reducing to the Wigner equation [17| without spin 



effects, which can be further reduced to the classical Vlasov equation. Here we start from the 
long wavelength limit (spatial scales much longer than the thermal de Broglie wavelength) 
of the full theory {3^, where all quantum mechanical effects are directly associated with the 
spin. A formal expression for the nonlinear low-frequency current is defined for a general 
geometry. This expression is then evaluated for two special cases with specified wave- vectors 
and polarization of the high-frequency fields. It is then found that although the classical 
terms gives the largest contribution to the nonlinear current for a plasma of moderate density 
and temperature, the spin terms can still contribute significantly to the generation of quasi- 
static magnetic fields. 

The organization of the paper is as follows. In Section II the calculation procedure 
starting from the classical Vlasov equation is outlined, and the nonlinear low-frequency 
current corresponding to this case is presented. In Section III the spin kinetic equation is 
introduced, and our main results starting from this model is derived. In Section IV we 
investigate the generation of quasi-static magnetic fields by the nonlinear current densities 
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that we have derived. In section V, our results summarized and discussed. Finally, the 
Appendix addresses the possible significance of quantum mechanical effects left out in the 
previous calculations. 



II. CLASSICAL KINETIC MODEL 



In order to outline our method of calculating the nonlinear current, we shall first consider 
the nonlinear mixing of two high-frequency waves (ci;i,ki) and {u2,^2) in a collisionless 
plasma. Considering only electron motion we describe the evolution of the perturbation / 
of the electron distribution function by means of the classical Vlasov equation 

dF„ 



f + v.V/+«(E 

ot m 



V X B 



-^(E + vxB).^ 
m av 



(1) 



where q/m is the electron charge to mass ratio and -Fo(v) is the unperturbed velocity dis- 
tribution function normalized such that uq = J Fq(Pv. For simplicity we assume that Fq is 
Maxwellian, i.e. Fq ~ exp(— f ^/2f^). Relativistic effects are neglected. For notational conve- 
nience we assume that U2 is negative. The nonlinear mixing between the two high-frequency 
waves thus yields a low-frequency response at (w, k) where u = Ui + UJ2 and k = ki + k2. 
Replacing the right hand side of equation ([T]) by (Ei + v x Bi) ■ {df2/d\) + (1^2) where 
(1^2) means interchange of indices 1 and 2, and /2 ~ — (zg/m)E2-(9-Fo/9v)/(a;2— k2-v), we 
then obtain the nonlinear and classical part of the generated low-frequency current density 
Jci as 



T — ^ 
Jci — 



iq 
m 



(cj — k ■ V 

which can be directly rewritten in the form 



El + V X Bi) ■ 1^ + (1 ^ 2) 
av 



dv 



(2) 



'cl 



q" 



k- V 



[Ei + v X Bi) ■E2 



El ■ VE2 ■ V ^ k2 



{UJ2 - k2 ■ V) 

(El + V X Bi)E2 ■ V 



;i^2) 



(3) 



Vf{uJ2 — k2 ■ v) ' {UJ2 — k2 ■ V 

Using the fact that the phase velocities of the high-frequency waves are much larger than 
the thermal velocity Vt we can adopt ki^2 ■ ^/^i,2 and uj/uji^2 as expansion parameters and 
obtain, up to order {uj/uji,2Y 



The first term of Eq. (jlj) (proportional to Ei ■ E2) is one order larger in the expansion 
made, but other terms are important as they are the ones contributing the generation of 
quasi-static magnetic fields. 



III. SPIN KINETIC MODEL 

Next we consider an extended distribution function /(r, v,s,t), where s is a spin vector 
with fixed length, describing the orientation of the electron spin. On semi-classical grounds 
it can be shown 29] that df/dt = 0, where d/dt is a total time derivative following the 
generalized particle orbit, i.e. including the spin evolution. Using the Heisenberg equation 
of motion to get ds/dt, a semi-classical model including the spin degree of freedom is found 



29|. With certain corrections, such a model can be derived as the long scale limit (scale 
lengths much longer than the thermal de Broglie wavelength) of a fully quantum mechanical 
treatment js^. The evolution equation in this regime is then found to be 



^(E + vxB) + ^V(s-B) + ^V(B-Vs) 

m m 2m 



■Vv/+^(sxB)-VJ = 



(5) 

where we here use the convention that |s| = 1. The difference from the semi-classical theory 
deduced from df /dt = is twofold. Firstly, the constant coefficients of the spin terms differ 
somewhat (this is compensated by different coefficients in the spin current given below, such 
that also the semi-classical model is energy conserving). Secondly, the term proportional to 
B ■ Vs has no semi-classical counterpart. This term can be viewed as part of the magnetic 
dipole force, which account for the fact that a quantum mechanical probability distribution 
of a single particle spin is always smeared out, as opposed to a classical spin (or magnetic 
dipole moment) which has a unique direction. The spin kinetic model is completed by 
Maxwell's equations with the current density 



,l.0.,n.^,jM. (6) 



where Jf and is the free part and magnetization (spin) part of the current density, 
respectively. Using spherical coordinates for the spin angles we have dQs = sinOsdOgdips- 

Next we follow the steps of section II, and write the nonlinear high-frequency source 
terms (due to (cji, ki) and {002, ^2)) on the right hand side. Since the zero order distribution 



function Fq is isotropic and independent of s we find 



dtf + v-Vf + 



^E+%V(s-B) 

m mn 



VvFo 



m 

2/ie 

h 



fE + V X Bl 



s X B) ■ V J. 



mh 



Vfs-Bl 



2m 



V(B-V« 



(7) 



Following the same procedure as in the classical case we consider a Maxwellian background 
distribution and use the notation 



k = ki + kg. 



(8) 
(9) 



Since the zero order distribution function is isotropic by assumption, the first order distri 
bution function can be approximated as 

-I (gEi,2 + ^V(s-Bi,2)) -^vFo 



and we get the free nonlinear current density 



(10) 



a; — k ■ V 



^(Ei + vxBi) + ^V(s-B) + |^V(B-V. 

m mn 2m 

2/ie 



h 



(s X Bi) ■ dj2 + (1 ^ 2) L 



(11) 



and the nonlinear magnetization current 

s 



iV X / d'^vdQ. 



a; — k • V 



m 

2/ie 



h 



El + V X Bi^ 
> X Bi) ■ dj2 + (1^2 



^V(s.B) + ^V(B.V. 

mn 2m 



■d.f2 
(12) 



A general geometry with arbitrary polarizations and wavevectors of the high-frequency 
waves leads to extremely lengthy algebra, and is not within the scope of the present pa- 
per. Thus as our first special case we consider parallel polarization and propagation. For 
definiteness we choose the geometry 



kl,2 = ki,2Z, 

Bi,2 = -Bi,2y, 

El,2 = -EloX. 



(13) 
(14) 
(15) 
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The high-frequency waves may be co-propagating or counter-propagating depending on the 
sign of ki^2- In this case subtracting the classical contribution Jd we obtain that the free part 
and the magnetization part of the nonlinear current due to the spin, denoted Jgp combines 
to 

Jsp = - q^hhB.B, [ dhdn, Fo 

m^f^ J uj — k v 



(16) 



{u2-k2Vz) V UJ2-k2Vz) {ui - kiV^) V "Vt UJi-kiV^ 

which for Vt ^ Wi / /ci , a;2/ A;2 reduces to 

1 A "^z k2Vz\ 1 A v1 kiV 



Jsp = -q^k^k2B^B2 / d'vdn 



m'^vf J o; — k ■ V 1002 \ UJ2 J ^1 \ Ui 

(17) 

Since k is along z, we note that Jsp in (fT7|) will be in the z-direction. This is the same 
direction as the classical contribution found from (jl]). Furthermore, for moderate tempera- 
tures and densities the classical contribution will be larger than that due to the spin, and 
hence the result given by ( ITTl) is merely a small correction, at least if we limit ourselves to 
parameters of temperature and density corresponding to laboratory and space plasmas. 

Next we modify our special case and consider the polarization of the waves to be per- 
pendicular. As it turns out, the spin contribution to the nonlinear current then vanishes, 
although there will be a significant magnetization in the z-direction. If we also modify the 
wavevectors slightly, this nonzero magnetization contributes to a nonlinear current density. 
This means that we consider the following geometry: 



El = Eix, 


(18) 


E2 = —E2y, 


(19) 


ki = kiz, 


(20) 


k2 = A;2Z + A/C2X, 


(21) 


Bi = 5iy, 


(22) 


B2 = 52X + AS2Z. 


(23) 



To limit the algebra, the deviation from parallel propagation is considered to be small, and 
we will thus only consider terms to first order in Ak2 or Ai?2- By Faradays law we note that 
AA;2//c2 = AB2/B2. When we calculate the nonlinear currents in this geometry it turns out 
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that the free part of the current vanishes, i.e. Jsp 
density reduces to 



' sp— m 



and the nonhnear spin current 



'sp 



mv?h 



UJ 



k- V 



f + 1 



k2 



CJl — ki ■ V CJ2 — k2 • V 



(24) 



which vanish if we let Ak2 — )■ 0. Taking the hmit of low temperature, i.e Vt ^ u/k, this 
further simplifies to 



'sp 



mh 



ki f k 



LOi \k2 



UJ2 \k: 



k 



B1B2 



(25) 



The geometry in this special case is of more interest, as we will find that the spin-contribution 
can give a larger contribution to the generation of quasi-static magnetic fields than the free 
current, also for relatively modest plasma temperatures and densities. This issue will be 
investigated further in section IV. 



A. A general isotropic background distribution 



The above results can easily be generalized to be valid for any background distribution 
that is a function of v'^, i.e. isotropic, which is of interest e.g. for a dense plasma when 
the thermodynamic equilibrium distribution is a Fermi-Dirac distribution rather than a 
Maxwellian. With this more general background distribution we get in the case of parallel 
polarization the nonlinear spin current 



Jsp =q^kik2BiB2 l^vdVL 



m 

X 



{UJ2 - k2V^) 



1-^ 

V? 



w — k ■ V 

^ k2Vz 



0J2-k2vJ (Ui-kiVz) 



1-^ + 



kiVz 



uji - kiv. 



d 



reducing to 

3sp = q^kik2B,B2 jdhdns 



cj — k ■ V 



002 



- 1-^ + 



k2Vz 



UJ2 



-11-^ , kr- 



(26) 

d 



(jJi J \ (9(f^) 
(27) 



for the limit Vt ^ uji/ki,U2/k2- 

In the case of orthogonal polarization in the same way we obtain 



'sp 



.32^ 
mh 



TTAk2 d^vBiBr 



(jj — k ■ V 



Wi — ki ■ V (jJ2 — k2 ■ V 



d 



(28) 



and in the low temperature limit this reduces to fl25|) . 



IV. MAGNETIC FIELD GENERATION 



To demonstrate the significance of the spin contribution in the nonhnear current density 
we write Ampere's law as 

ik X B{uj, k) =/io (Ji(a;, k) + Jju, k)) - ^E(u;, k). (29) 

Next we write the linear current density as Ji(a;, k) =(t{u, k)E(a;, k), where a{uj, k) is the lin- 
ear conductivity. Combining (l29l) with Faraday's law, we immediately find the low frequency 
magnetic field generated by the nonlinear current as 

z^pk X z/ipk X (J^i + Jsp) 

^) = D:jij:^ = DUu^,k) 

where Dcm{^, k) = — k'^c^ — icua^co, k) — k'^c^ — cUp. Eq. fl30|) thus contains both the 
classical contribution from (jl]) as well as the spin contribution from flMl) . To shed further 
light on this expression we evaluate ( 130|) in the low temperature regime kvth <^ i^, in which 
case the expression simplifies to 

iHok BiB2Ak2 ( q^v'fikj + k^) IG/i^vr/ca [^i ,\ ^2 f k 

B(w,kj = -— — - — tt'^o- 



Dem{^,k) m 002 \ mkik2Ui2 ^ 



(31) 

Thus we see that the spin contribution to the magnetic field generation (proportional to 
fJ'l/h) dominates in the regime 

h'k,k2 > ^ (32) 

whereas the classical contribution dominates otherwise. A case of experimental interest 
clearly includes two high-frequency sources. However, we can note that our results also 
are of relevance for a single source, where (ciJi,ki) and (ci;2,k2) represents different spectral 
components of a focused pulse. In this case typically |k| ^ |ki^2|- For the condition ( l32l) it 
does not matter whether two pulses or a single source is used. Furthermore, for current laser 
plasma experiments with lasers in the optical regime, it is clear from ( !32|) that the classical 
contribution will dominate. However, in case laser-plasma experiments with an X-Fel source 
such as that being built in DESY is made js^j, |ki,2| ~ 6 x 10^m~^, and the quantum spin 
effects can be of importance for magnetic field generation for plasma temperatures T < 2 
xlO^K.This condition have been derived for the case of a specified geometry, and it should 
be noted that the results may differ in case polarizations and/or the directions of wave 
vectors are changed. 



V. SUMMARY AND DISCUSSION 



In the present paper we have studied the nonhnear current density generated by high- 
frequency waves in a plasma, with a focus on the contribution emanating from the electron 
spin. The largest part of the current density is usually associated with the classical pondero- 
motive effect. However, it is found that although the largest part of the nonlinear current 
in a moderate density, moderate temperature plasma is due to the classical terms, the spin 
may still give a significant contribution to the magnetic field generation mechanism. For 
the geometry considered here, the condition needed for spin effects to be important require 
short-wavelength sources, of the order of the x-ray regime. Besides the quantum effects due 
to spin considered here, there is also particle dispersive quantum effects. Although a thor- 
ough consideration of such effects is still to be made, our calculations outlines here indicate 
that particle dispersive terms may be of comparable importance for the nonlinear current 
density. Thus there still remains much research in this area to be made. 
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Appendix: Particle dispersive efTects 

The quantum effects due to spin that we have considered here should be compared to 
the ones due to particle dispersive effects. Such effects can be described by the Wigner 
function, that reduces to the classical distribution function whenever the thermal de Broglie 
wavelength is small compared to the scale lengths of the problem. The quantum corrections 
due to this in the kinetic equation scale as h"^, whereas e.g. the magnetic dipole force due to 
the spin scale as h. However, this does not necessarily mean that the lowest order quantum 
corrections is due to the spin, since for a spin independent zero order distribution function, 
the spin term at one place (e.g. the magnetic dipole force) always need to be combined with 
another spin effect (e.g. the magnetization current) to produce a non-vanishing nonlinear 
current. Thus both the spin effects and particle dispersive effects produce quantum correc- 
tions that are proportional to h"^ to lowest order. While our main focus of the present paper 
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are the corrections due to the spin, we will here briefly outline how to obtain the quantum 



correction due to particle dispersive effects. The general e vo 
spin and particle dispersive effects was derived in Ref. 



ution equation containing both 



30| . Here we take that evolution 



equation, drop the spin effects considered above, and keep the particle dispersive effects in 
a weak quantum expansion (with the thermal de Broglie wavelength over the characteristic 
scale length as expansion parameter), where only terms up to are kept. The governing 
equation then reads 

7^/ + v-V,/ + ^(E + vxB)- V,/ 

ot m 



2Am? I m 



— (E + V X B 



— B X V 

m 



V (yx-^v)] ■ (^B X V, + V, 

(A.l) 



/ 



(A.2) 



Deflning 

f = Fo 

where Fq is the background distribution and / is the perturbed distribution function, we 
can separate ([7]) in linear and nonlinear terms: 



^/ + V ■ V,/ + ^ (E + V X B) ■ V.Fo 
ot m 



E + V X B ■ V, 



^(E + vxB)- V./ 



m 
2 



^|^(E + vxB 



24^3 -"(^-^^ 



-B X V. f^. -^v)]- f-B X V.Fo + V./) 1 
m V /J \m J J 



(A.3) 



Assuming two waves as in previous calculations, and a Maxwellian background distribution 
we have 



ig(Ei,2 + V X Bi^a) ■ V^, 



24m2 



/i,2 ~ -, S; ^ — • A.4 

where we separated the perturbed distribution function into a linear part /i 2 and a nonlinear 
part jni- The rest of the calculations can be performed as in section III, although the algebra 
gets extremely complicated in general. A thorough treatment of particle dispersive effects is 
beyond the scope of the present paper, but we will nevertheless point out two conclusions. 
Firstly, that the particle dispersive effects can be comparable in magnitude to the spin 
contributions. Secondly, although these two quantum effects may be comparable, they do 
not typically cancel, as the spin current for parallel propagation parallel polarization scales 
as oc W'k\k2BxB2 (see eq. (|T7|) ). whereas the corresponding scaling can be shown to be 



/ 
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oc h'^BiB2/kik2 for particle dispersive effects. Thus we conclude that the contribution of 
particle dispersive effects to the nonlinear current density is worthy of consideration in its 
own right, but we should not expect such contributions to cancel those due to the spin. 
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